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Abstract 

\q , We discuss Staruszkiewicz 's nonlinear modification of the Schrodinger equation. It is pointed 

out that the expression for the energy functional for this modification is not unique as the field- 
theoretical definition of energy does not coincide with the quantum-mechanical one. As a result, 
this modification can be formulated in three different ways depending on which physically relevant 
properties one aims to preserve. Some nonstationary one-dimensional solutions for suitably 
chosen potentials, including a KdV soliton, are presented, and the question of finding stationary 
solutions is also discussed. The analysis of physical and mathematical features of the modification 
leads to the conclusion that the Staruszkiewicz modification is a very subtle modification of the 

£^ ■ fundamental equation of quantum mechanics. 
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1 Introduction 



Some time ago Staruszkiewicz |T] put forward what seemed to be at the time unique^ for three dimen- 
sions a way of modifying the Schrodinger equation by adding to its action density a Galilean invariant 
term Aw = l(AS) 2 /2, historically motivated by his theory of the free electromagnetic phase |2|, |3|. 
Here S stands for the phase of the wave function \l/ = Rexp(iS), originally identified with the electro- 
magnetic one. The coupling constant 7 was found by some reasoning involving electromagnetic theory 
to be equal 1/167T 2 in natural units 0. This constant is dimensionless in the system of natural units 
in three dimensions, a property that distinguishes the Staruszkiewicz modification from any other 
proposed so far. However, neither the identification mentioned nor any other arguments based on 
electromagnetic theory are really necessary as the modification in question can simply be postulated 
on the grounds of physical consistency.^ In such a case 7 is an arbitrary constant. It is this approach 
that we adopt for our paper. It should be stressed that this is a purely field-theoretical approach. As 
of now the equations proposed by Staruszkiewicz have not been derived within any phenomenological 
scheme. Unfortunately, this is rather disadvantageous for the physical analysis of the modification as 
any insight one could derive from such a scheme cannot be gained here. 

The purpose of this paper is to extend the original brief report ]]J by providing a more elaborate 
analysis of some aspects of this interesting modification together with some solutions to it. To this 
end, in the next section, we review main points of the modification largely along the lines of and 
discuss them in some detail. The following section contains a more detailed study of the modification, 
including a comment on the variational derivation of the equations of motion, a discussion of the 
problem of separability of noncorrelated subsystems, remarks on its semiclassical approximation and 
the Ehrenfest relations, as well as other comments on its characteristic properties. This section also 
addresses the problem of the energy definition, which for nonhomogeneous nonlinear modifications 
as the one in question is not unique J7|. Depending on how one chooses to approach this issue one 
can have two different formulations of the modification. Yet another formulation can be adopted to 
cirmcumvent the problem of nonseparability in the approach of Bialynicki-Birula and Mycielski [|10] . In 
another section, concerned with solutions, guided by earlier observations, we find some characteristic 
and physically important solution, a coherent state. It is in this section that we also demonstrate 
how to reduce the continuity equation to the form of the Korteweg-de Vries equation and the entire 
Schrodinger equation to a soluble form for a suitably chosen but time-dependent potential. Here we 
also address the problem of additional, beyond those that satisfy the linear Schrodinger equation, 
stationary solutions for the one-dimensional case. Other solutions can be found when the equations of 
motion of the modification are cast into a dynamical system, a stationary point of which turns out to 
be a plane wave. The coherent state solution and the Gaussian wave packet are also identified within 
this framework. Our findings and observations are summarized in the last section, where a broader 
motivation for this study is also presented. 



x It was demonstrated in g] that the Staruszkiewicz modification can be extended by allowing more Galilean invariant 
terms of the same property in three dimensions. 

2 In fact, as shown in |6) any fully electromagnetic extension of the Schrodinger equation based on the identification 
of these phases runs into conflict with the Galilean invariance of the modified equation. 



2 



2 Main Features of the Modification 



To begin with, let us recall after Q that with the term 7(AS') 2 /2 implemented the modified action 
for the Schrodinger equation 
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The other equation, that remains unchanged in this Schrodinger-Madelung representation, reads 
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The first of these equations should be viewed as the continuity equation for the probability density 
p = R 2 and the current 
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The solutions to (2) and (3) are supposed to be normalized in the sense / d 3 x R 2 = 1, unless they 
correspond to the continuous part of the energy spectrum, in which case the normalization by imposing 
the periodicity condition, or the normalization to the delta function is to be applied. Along with this 
we need to consider the energy functional which in the original formulation of Staruszkiewicz M was 
defined as 
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This is a field-theoretical definition of energy. One arrives at it from a given Lagrangian density L, by 
identifying the energy density of a field configuration with the time-time component of the canonical 
energy-momentum tensor 
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Here L is assumed to depend on a set of fields fa and their first and second order derivatives. The 
total energy Eft is then given by a space integral over Tq , one example of which is formula (5). The 
energy-momentum tensor satisfies the conservation law, 
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resulting in the energy being a constant of motion unless the parameters of Tq (like the potential V) 
depend explicitly on time. In some cases though even such time- dependence does not prevent Eft 
from being a constant of motion. 

The Staruszkiewicz modification can be put in the form that involves the entire wave function 
\1/ and not only its parts R and S. By observing that AS = — iAhx (^ r */^ r ) /2, one can write the 
Lagrangian density for the modification as 
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The factor 7/8 in the last expression was chosen so as to reproduce the continuity equation in the 
form of (2). By altering the density in question so that the modification term is homogeneous in the 
wave function \l/ similarly as the rest of the Lagrangian one arrives at yet another modification ||. 
Comparing these two almost identical modifications can provide a good way to learn about the impact 
of nonhomogeneity. From (8) one obtains now the modified Schrodinger equation for \1/ 
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As seen from this equation, the Hamiltonian for the Staruszkiewicz modification, Hsm, can also be 
couched in the hydro dynamic form as 
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where H^se is the Hamiltonian for the linear Schrodinger equation. Knowing that in quantum theory 
the energy of a quantum system is defined as the expectation value of the Hamiltonian H, that is, 
E QM =< yff\H\V >= J d 3 x^f*H^, one finds that Q 
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Now, even if these two energy functionals are equal for most solutions to the linear Schrodinger 
equation that are also solutions to this modification, they do drastically differ for ordinary Gaussian 
wave packets for which AS = g{t). Noting that Eq M is not the same as E FT , one is faced with a 
challenge of choosing the right form of the expression for energy. We will discuss this issue in the next 
section. 

It is reasonable to demand that E < 00, which besides the normalization condition constitutes 
another constraint on the solutions to (2) and (3). In fact, as just pointed out after |]] , this constraint 
cannot always be satisfied, even for solutions as common to the Schrodinger equation as the Gaussian 
wave packets if Eft is adopted for the definition of energy. 

To conclude this section, let us observe that without the modification term equations (2-3) con- 
stitute a linear equation, the ordinary Schrodinger equation, in a disguised nonlinear form. One can 
always transform them back to the linear equation by employing S(^f, \&*) = -iln($*/^) /2 and 
R 2 (^/, = \E f \E'*. On the other hand, the term A 2 5* is linear on its own. Therefore, paraphrasing J. 
A. Wheeler, one could say that the Staruszkiewicz modification contains "nonlinearity without non- 
linearity" . It is, in some sense, a minimal way of introducing nonlinearity and perhaps unique to this 
modification in which the linear equation in the disguised nonlinear form is supplemented by another 
linear term. The latter however becomes nonlinear in the Schrodinger representation where A 2 ^ turns 
into — iAln (\I/*/^) /2- This incidentally indicates that this nonlinearity may not be linearizable, i.e., 
it is not reducible to the linear Schrodinger equation by means of a suitable transformation. 

As noted, the discussed modification does not introduce any new dimensional constants in the 
system of natural units in three dimensions. In an arbitrary system of units the dimension of the new 



3 The continuity equation implies that for all those situations where R 2 is normalizable and the gradient of the 
phase decays fast enough at infinity, the imaginary term in (11) does not contribute at all by the virtue of probability 
conservation. 
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coupling constant 7 is the same as that of the product he or the square of the electric charge. The 
coupling constant 7 can thus be represented as 7 = 'y'hc, where 7' is a dimensionless parameter in an 
arbitrary system of units and c is the speed of light. In a way, if we are to paraphrase J. A. Wheeler 
once more, this product can be viewed as introducing the self-interacting electric charge without 
introducing any charge at all. However, the combinations Gm 2 and V Ghcm, where G is Newton's 
constant, have the same dimension as he. This suggests that 7 could be of the order of either of the 
discussed combinations, meaning that 7' is of the order of 1, although diametrically different values of 
this parameter cannot be excluded by these purely dimensional considerations. In the case 7 is of the 
order of Gm 2 or V Ghcm, one can suspect that the Staruszkiewicz term (AS) 2 is of a gravitational 
origin. 

However, the dimensionless nature of the coupling constants of the Staruszkiewicz modification is 
by no means universal for it is only characteristic of three dimensions in the natural system of units. 
One can also argue whether the choice of the system of units in which c = 1 is physically justifiable in 
a modification of the nonrelativistic equation and whose nonrelativistic status is sealed by its Galilean 
invariance. We believe that similarly as the original electromagnetic underpinnings of the modification 
cannot be maintained for this would violate the ?7(l)-gauge invariance and the Galilean invariance 
of the scheme 0, so it is more consistent to treat the Staruszkiewicz model without referring to any 
special system of units, particularly that this is not invariant. Let us note at this point that there 
exist quantities that are dimensionless in any system of units. Such quantities, to name as an example 
the fine structure constant, can be defined as ratios of other quantities having the same dimensions. 
Yet, the coupling constant in the Staruszkiewicz proposal is not of this nature. These shortcomings of 
the original motivation of the Staruszkiewicz proposal can easily be avoided. To this end, we suggest 
that this proposal be reformulated as the simplest nonhomogeneous phase modification that derives 
from a local Galilean invariant Lagrangian. The assumption that it is the Lagrangian that should be 
Galilean invariant and not only the equations of motion together with the assumption of simplicity 
leads to a unique modification with a dimensionless coupling constant in the natural system of units 
in three space dimensions. It should be noted that in order for the simplicity in question to become 
apparent, the modification needs to be put in the framework of hydrodynamic representation of the 
Schrodinger equation. Reformulated in this way, the Staruszkiewicz modification can truly stand on 
its own merits. 



3 Further Analysis of the Modification 



It may be instructive to see how the equations of motion for R and S are derived for the modification 
under study. The nontrivial aspect of this derivation consists in the fact that the action functional (1) 
from which these equations are obtained contains derivatives of the second order and because of this 
it yields a new boundary term, coming from the (AS 1 ) 2 part of the action. Varying (1) with respect 
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where D stands for the domain of integration, being the four-dimensional space-time, and dD for its 
boundary which is spatial infinity in the limit \t\ — > oo. The second integral in the last expression 
is the new boundary term. To ensure that SS mo d = one needs to assume either that 5(VS) = 
on the boundary or that AS vanishes therefj The latter choice seems to be more attractive, as one 
would expect it to promote solutions with finite energy. Unfortunately, even if AS falls off sufficiently 
rapidly for large \x\ or large t, the energy functional (5) or (11) can still have singularity at the origin. 
Therefore, in general it seems to be unavoidable to formally impose 5(VS)\dD = 0, which in fact may 
constitute a constraint on the class of acceptable functions for S. In any case, in order to arrive at the 
Staruszkiewicz modification it is necessary to ensure that the discussed integral does not contribute 
to the variation in question. Of course, if the equations of motion for the modification are simply 
postulated and not derived from a Lagrangian density, the above analysis does not apply. 

The fact that the Staruszkiewicz modification is not homogeneous causes it to be nonseparable in 
the sense first discussed by Bialynicki-Birula and Mycielski [TU|, also referred to as the weak separa- 
bility. This approach is valid only for noncorrelated subsystems. That this is so was demonstrated 
for its extended version in |5|]. The issue of separability is certainly physically important. The lack of 
separability means that even in the absence of interactions the motion of one wave packet can affect 
the other in a system consisting of these two packets. This question has recently been reformulated 
by Czachor |l l]| in a novel "effective" way which, instead of pure states, uses density matrices as 
basic object subjects to a nonlinear quantum evolution. It has been shown that a much larger class 
of nonlinearities is allowed than in the Bialynicki-Birula and Mycielski approach. In particular, all 
nonlinearities of the form F(\ip(x)\), as for instance the cubic nonlinear Schrodinger equation, are ac- 
ceptable. As opposed to the weak separability that assumes that the total wave function of the system 
is factorizable to reflect the fact that its subsystems are noncorrelated, the wave functions in Czachor's 
approach do not have to obey this condition. For this reason, the approach in question belongs to 



the category of strong separability. As asserted by Czachor |TTJ, the Staruszkiewicz modification is 
separable in this framework. 

As noted in the preceding section, the energy definition for the Staruszkiewicz modification exhibits 
a curious ambiguity that is absent in the energy formulation for the linear Schrodinger equation. Since 
quantum mechanics is a probabilistic theory, it is certainly reasonable to adhere to the probabilistic 
interpretation of energy as a quantity that is defined only in terms of averages. Therefore, the correct 
form is given by the expression for Eqm which as the mean value of the Hamiltonian represents the 
energy according to the standard interpretation of quantum theory. In keeping with this approach, 
one dismisses the integral of motion derived from the conservation of energy-momentum tensor as 
a viable candidate for energy in the quantum-mechanical framework, although its predictive power 
can still be useful for systems of finite Eft- As a consequence, since Eft is not really an observable 
in quantum theory unless it matches the expectation value of the Hamiltonian, one should be not 
concerned if this quantity is infinite, even though any infinity in a physical theory is rather disturbing. 
A less radical alternative approach would consist in reconciling these two different objects. 

Let us now argue that compared to Eft, Eqm seems to be more acceptable on the grounds of 
physical consistency. In the Staruszkiewicz modification the phase is decoupled from the amplitude 
of the wave function for, as seen from the equations of motion, it remains a physically determined 
quantity, having a life of its own even if R = 0! In linear quantum mechanics, except for the nodal 



4 The other boundary terms in I2 are handled in the standard way by assuming that variations of both R and S 
vanish at spatial infinity for \t\ — > 00. 
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points that constitute a discrete set of points, a situation like that does not take place. In fact, it 
seems to be meaningless to talk about the phase when the amplitude vanishes in a finite region of 
space and so do the equations of motions leaving no room to determine the evolution in this region. 
However, in the Staruszkiewicz modification the phase can, in principle, be determined through the 
equations of motion even if it is not accompanied by the amplitude. To see what consequences this 
can possibly entail let us consider a one dimensional wave packet moving freely in space until it 
encounters a totally impenetrable wall, an infinitely large potential barrier. Now, according to the 
standard interpretation of quantum mechanics, the amplitude of the packet has to vanish beyond the 
wall, but this does not apply to the phase which can be nonzero there. In linear quantum mechanics, 
the energy associated with the phase beyond such an infinite wall is zero as long as its gradient is finite 
which is to be expected as the probability current j = R 2 VS should vanish there, too. This is not so in 
the Staruszkiewicz modification unless one requires that the energy contribution due to the nonlinear 
part is zero in the region beyond the wall. The vanishing of the probability current does not entail this 
if one identifies the energy with Eft- Therefore, one is faced with a transmission of energy through 
an infinitely large barrier that may not be accompanied by a flow of the probability current. It is clear 
that this situation is not physically sound, but fortunately it can be amended if one defines energy as 
Eqm-i which seems to be the first argument in favor of this definition of energy. Now, assuming that 
the current is a continuous function and vanishes on the wall implies that it vanishes also beyond it as 
A 2 5* is zero everywhere beyond the wall by the virtue of the equations of motion or just the continuity 
equation. In other words, A 2 S = implies that the current is constant beyond the wall, but owing 
to its continuity it can only be zero there. Moreover, the energy associated with the phase beyond 
the wall is zero as is the energy of a "quantum mechanical system" in the entirely free space in the 
absence of the amplitude if identified with Eqm- Unless the phase alone is really proven to be endowed 
with energy, this observation once again supports Eqm as a more physically reasonable expression 
for the energy of a quantum-mechanical system. Still, even if no energy can be associated with the 
phase in the limit of vanishing amplitude, it is not out of the question that it can be detected by some 
diffraction or interference phenomena similar to the Aharonov-Bohm effect [|13f . To summarize our 
reasoning, electing Eqm over Eft stems from a very simple yet physically respectable requirement 
that no information, energy in particular, is allowed to be transmitted through an infinitely large 
potential barrier. 

The energy functional Eqm contains an imaginary component. Since the energy is supposed to be 
a real quantity one might want to require that this part does not contribute to the total energy, which 
imposes a constraint on physically acceptable states in this particular nonlinear model of quantum 
mechanics. That this constraint is not necessarily very restrictive has already been demonstrated for 
the Gaussian wave packets. Since for such packets Eft is infinite, they would have to be excluded from 
legitimate solutions to the modification if this energy definition were employed. Yet, they are perfectly 
fine on the energetic grounds if one uses the quantum-mechanical definition of energy. Moreover, the 
continuity equation coupled with the assumption of probability conservation implies the vanishing of 
this term for a large class of physically plausible situations. The proposed condition of real energy 
is equivalent to the selection of observables in nonlinear quantum mechanics for which contributions 
of nonlinear parts vanish on normalized states [12]. On the other hand, complex energy is not in 
the least foreign to linear quantum mechanics and therefore it is certainly justifiable to entertain the 
possibility of retaining it in the modification under discussion. In linear quantum mechanics, complex 
energy can appear as a result of implementing complex or "optical" potentials usually invoked to 
phenomenologically model absorption in scatttering processes. The complex potentials have also been 
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used to describe decoherence |TJ]. Since the Staruszkiewicz modification introduces complex potential- 
like terms in a natural manner, one finds it tempting to allow states with complex energy if only in 
the hope that they can mediate in the process of decoherence. Such states, however, would have to be 
rather exotic because, as noted, in the majority of physically reasonable situations one expects the 
integral over the imaginary term in question to vanish. 

As we see, the properties of observables depend on the space of states in which they are defined, 
which also applies to linear quantum mechanics, where, for instance, the self-adjointness of an operator 
depends on its domain. In line with this approach, one can attempt to reconcile Eqm with Eft by 
choosing a domain in which they are equal for each function in the domain. This alternative seems 
to be coming as close as possible to the original Staruszkiewicz formulation by retaining Eft as a 
physically relevant object, but it is probably more restrictive than the formulation just proposed for 
it requires not only that Eft be finite and Eqm be real, but, in addition, that these two quantities 
be equal. In what follows, for the sake of exploring all possibilities, we will consider both of these 
formulations. In fact, none of them seems to be decidedly more convincing than the other. We will 
refer to the alternative in question as maximal to contrast it with a moderate formulation that is 
concerned only with the expectation value of a Hamiltonian treating Eft as useful but not a crucial 
quantity in the quantum-mechanical scheme it adheres to. As we will see, all the nontrivial solutions 
that we have found belong to the moderate formulation, although one of them is shared by the 
original formulation as well. This lack of uniqueness in the definition of energy clearly marks the 
difference between linear quantum mechanics and its nonlinear models where uniqueness seems to be 
an exception rather than a rule. 

It is possible to adopt yet another approach that takes care of the weak nonseparability problem 
and the problem of uniqueness of energy at the same time. This is accomplished by considering 
as legitimate only those solutions to the modified Schrodinger equation that satisfy the equation 
A 2 ,S = 0. Even if the strong separability seems to offer a more all-encompassing approach, it is 
the weak separability that is a more stringent condition and for this reason one can use it as a 
criterion for the formulation discussed. We will call this formulation minimal for the domain of its 
solutions is expected to be the smallest of all the formulations considered. It overlaps with that of the 
maximal formulation in the case AS = and is entirely contained within the domain of the moderate 
formulation. Despite its name, this approach ensures the maximum of physically desirable properties, 
including, as we will soon see, the standard classical limit of quantum theory. 

Let us now comment on the semiclassical approximation for the modification. Formally, one 
performs this approximation by rescaling S to S/h, followed by expanding S in a series S = So + 
%S\ + /i 2 S l 2+...One obtains in this way an expansion of the Schrodinger equation in the powers of 
the Planck constant. The lowest order approximation term to the continuity equation of the linear 
Schrodinger equation is proportional to %. This, however, is not the case for its modified counterpart 
whose expansion in the powers of h starts with h~ l . The next leading term is proportional to % . 
As a result, we obtain that A 2 So = A 2 Si = 0. In practice, this means that, in a very good first 
approximation, the solutions to (2) can be sought for among wave functions whose phase S satisfies 
the condition A 2 S* = 0. One could hope to find in this way even exact solutions to (2) and (3), with, 
perhaps, some additional but physically well justified potential V. As we will see in the next section, 
this hope turns out to be justified. 

The classical limit of this modification in the sense of the Ehrenfest theorem may not always exist 
since the standard Ehrenfest theorem of linear quantum mechanics is altered by nonlinear corrections. 
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As worked out in [B], the general form of the modified Ehrenfest relations is 



mj t (r) = (p) + F u (12) 
|(P) = -(W) + F 2 , (13) 



where 

2m 



F 1 = ^J d 3 xfR 2 Hj, (14) 
d 3 xR 2 (2HjVS - VH R ) . (15) 



Hr and Hi represent the real and imaginary parts of the nonlinear part of the Hamiltonian of the 
modified Schrodinger equation, respectively. In the derivation of the last formula it was assumed 
that / d 3 xV (R 2 Hf) = and because of that this term was discarded. Since Hj = r yA 2 S/2R 2 , 
this is tantamount to assuming that IxpA 2 ^ vanishes at infinity. This, in turn, is implied by the 
condition that the energy Eqm be finite, and it is satisfied for all those configurations that are 
square-integrable as can be seen from the continuity equation given the probability is conserved. Had 
/ d 3 xV (R 2 Hj) been nonzero, it would introduce an imaginary component to F2. Now, since Hr = 
and Hj = , yA 2 S/2R 2 , the Ehrenfest relations for the modification in question read 

m l (r) = (p) + ^ / d 3 xfA 2 S, (16) 
at ' n J 

j f (p) = - (W) + 7 J d 3 xA 2 SVS. (17) 

These relations are Galilean covariant and their nonlinear contributions are also Galilean invariant for 
the wave functions that are normalizable and do not violate the probability conservation. 

We see that, in general, the nonlinear corrections to the Ehrenfest relations do not vanish, lead- 
ing to a different classical limit than in linear theory. This feature of the modification is shared by 
other nonlinear generalizations of the Schrodinger equation, the only notable exception from this rule 
being the Bialynicki-Birula and Mycielski modification. It suggests that the Staruszkiewicz modifica- 
tion is not linearizable and as such may describe phenomena that cannot be captured by the linear 
Schrodinger equation. On the other hand, perhaps due to its simplicity, the relevance of the Ehrenfest 
theorem tends to be overestimated. As argued in ||15| , this theorem is neither sufficient nor necessary 
to characterize the classical regime of quantum theory. 

However, one can also entertain the possibility that the classical limit of quantum world does 
not always exist in the form suggested by classical theory. This might offer an explanation of some 
anomalous or otherwise unexplained phenomena some of which owing to their strangeness remain 
outside the realm of accepted exact science. As a matter of fact, since even the quantization procedure 
is sometimes ambiguous, a unique classical limit of quantum theory is something that is enforced by 
rather than naturally emerges from it. What is remarkable in the modification discussed, it is that 
the departure from the standard classical limit is very subtle and, at the same time, pretty elusive 
due to the very elusive nature of the phase. The departure in question, unlike, for instance, in some 
other modification that possesses the same property [0], is also entirely defined by the phase.Q 

5 Other modifications, as for instance the Doebner-Goldin [Q modification, can also have a non-standard classical 
limit. However, since the Doebner-Goldin modification describes only irreversible and dissipative quantum systems it 
cannot be treated on the same footing as the Staruszkiewicz modification. Moreover, it is also less unique for it involves 
as many as 6 dimensional parameters. 
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Finally, let us make two more comments utilizing yet another rather straightforward modification 
of the Schrodinger equation, a relativistically modified Schrodinger equation which follows when the 
expression for the energy of a single particle in terms of its momentum is taken in the leading relativistic 
approximation 

tf=^---4r, (18) 

and the first quantization is performed. The Schrodinger equation for this modification, 

•H-^-SH?*)-* <19) 

when put in the Madelung representation involves additional contributions to the continuity equa- 
tion in both S and R. It is one of the reasons that this particular modification can be considered 
more general than the Staruszkiewicz modification. However, despite its generality, the relativistic 
modification does not contain Staruszkiewicz's modification as a special case. This can be seen by 
comparing the dimensions of the coupling constants of these modifications. Moreover, unlike the 
Staruszkiewicz modification, the relativistic modification is homogeneous of degree one in the wave 
function and is linear. The other comment consists in noting that the Gaussian wave packets, which 
are not solutions to the Staruszkiewicz modification for energetic reasons^], are also not solutions to 
(19), which indicates that it is rather difficult to preserve wave packet solutions in other extensions of 
this equation as well. Therefore, the modification under discussion in its original formulation should 
not necessarily be viewed as pathological in this respect. However, what is really unacceptable here 
is that Eft 7^ Eqm- This causes the wave packets to be excluded from the maximal formulation of 
the Staruszkiewicz modification. 



4 Solutions 



It is not easy to find solutions to the modification discussed here. The main obstacle to this stems 
from the nonhomogeneity of (2). Even if such solutions are found they may not necessarily yield 
to a clear physical interpretation. Therefore, it should not come as a surprise that the physically 
best understood solution we discovered so far represents a case when the nonhomogeneity does not 
interfere. 

As suggested in the previous section, one should look for the solutions to the Staruszkiewicz 
modification among the phenomena for which the square of the Laplacian of the phase S vanishes. 
We do not intend to present here a complete solution to the equation A 2 S = together with a 
discussion of its applicability to this modification, reserving it for another publication |J . Instead, 
we prefer to show the simplest possible and physically important solution to the equation AS = 
that, obviously, satisfies the previous one and ensures the finite unique energy in each of the discussed 
formulations of this modification. Such a solution is provided by a coherent state ||. It is easy to 
check that indeed the coherent state described by 
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6 As found in 
wave packets. 



1 , due to a slow convergence of their phase in the infinity, the energy Eft is not finite for the Gaussian 
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and 



S - [ — - ^ sin 2 (out - 8) + ^ a[X sin (out - 8) ) (21) 
2 x / 



represents a solution to equations (2) and (3) in one dimension in the potential of a simple harmonic 
oscillator V = mcu 2 x 2 /2. Here xq = (mu;) -1 / 2 , while a and 5 are arbitrary constants, complex and 
real, respectively. We also put % — 1. The three dimensional case can be worked out as a product of 
one-dimensional amplitudes with ou and the phase 5 appropriately redefined. 

Another physically interesting solution to (2) in one dimension can be found for a general negative 
7 = — 17|. In this case, (2) becomes 

dR 2 h 2 

h— + - (R 2 S X ) x + | 7 |^ = 0. (22) 
ot m v ' 

We will now show that this equation can be reduced to one of the standard forms of the celebrated 
Korteweg-de Vries equation 

u t + 3(u 2 ) x + u xxx = 0. (23) 

To this end, let us make identifications S x = R 2 = du(x), where d is a free constant whose dimension 
is to be worked out. Because the dimensions of S x and R 2 are equal, the dimension of t has to 
be meter 3 . This makes it impossible to maintain the natural system of units if we insist that 7 be 
dimensionless to preserve the fundamental idea of the modification in its original formulation. We can 
assume that either h = 1 or c = 1 but not both. Such a compromise is unavoidable in one dimension. 
We will assume that h = 1 which implies that the speed of light has the dimension of meter -2 and 
the dimension of m has to be that of meter. Now, it is certainly in the spirit of this modification to 
make the coefficient of (u 2 ) x dimensionless and, in fact, it has to be so for (22) to reduce to the form 
of (23). Therefore, d/m\ , y\ = 3 which does not impose any constraints on I7I because of the freedom 
of choosing d, but fixes the dimension of d to be that of meter. Moreover, t gets rescaled to t' = l^t. 
The most physically attractive solution to the Korteweg-de Vries equation is known to be a soliton of 
the form 

"■^ k) = ^i- (24) 

where 8 = + t, £ = x — |7|fc 2 £, and r is an arbitrary phase. The dimension of k 2 is the same 
as the dimension of u(x), meter -2 . This solution is normalizable and the normalization condition, 
/ dxR 2 (x) = dj dxu(x) = 1, implies that dk = 1/2. Consequently, since d/m\j\ = 3, the "speed" of 
the soliton is v s = |7|/4<i 2 = l/36|7|m 2 . From (3), we find out that the entire Schrodinger-Madelung 
system of equations can support the KdV soliton only in the presence of the potential 



V(9; k) ' 



8m cosh 4 % 



f) f) 
k 2 + AEm) cosh 4 - + 2A; 2 (2md\~i\k 2 - l) cosh 2 - - d 2 k 4 



(25) 



2 

where E is a constant of the same dimension as energy, meter -3 . In this potential Eft is finite (and 
also constant!), while Eqm is real and constant as well, however they are not equal. It is not clear how 
this potential can be physically realized. The solution under consideration, similarly as the Gaussian 
and coherent state wave packets, does not alter the standard Ehrenfest relations. 

The solutions discussed above were non-stationary. As observed in ||, the Staruszkiewicz modifica- 
tion, as opposed to its extension || and some other notable nonlinear modifications of the Schrodinger 
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equation |TI| [L8| , does not affect stationary states of quantum mechanical systems determined by the 
Schrodinger equation. This means that wave functions describing these states are also solutions to 
the modified equation with unchanged energy levels. It is this particular property that makes the 
Staruszkiewicz modification stand out. As a function of time and position, the phase of these states is 
given by S = —Et/h+cr(x), but its position dependence encoded in a{x) is physically inconsequential. 
Usually, a(x) is at most a linear function of spatial coordinates. However, it is conceivable that there 
exist other stationary states for which the phase is also a physically nontrivial function of position in 
the sense that, for instance, the phase affects the energy of the system. We will now show how one 
can approach the problem of finding such solutions in the one-dimensional case. 

To find stationary solutions to the Staruszkiewicz modification one starts with the condition 

dR 2 , , 

0. (26) 



dt 



This condition when applied to (2) and (3) together with the assumption of time-independent potential 
implies that dS/dt = const = —E/h, where E is to be identified with the energy of a quantum- 
mechanical system. For one-dimensional stationary problems, assuming h — 1, equations (2) and (3) 
reduce to 

u" + aR 2 u + c = 0, (27) 
R" + 2m (E - V - bu 2 ) R = 0, (28) 

where u = S x , a = — 1/7771, b = l/2m, and c is a free constant that can be zero. 

One can seek solutions to this system in two different ways. One is by selecting a normalizable R 
representing a solution to the linear Schrodinger equation in some potential V. Treating this as an 
Ansatz one attempts to solve the first of the above equations for u. This method will usually lead to a 
complicated equation for u that may not be solvable in an analytical manner even for a simple R. Once 
u is determined, the elected R can be interpreted as the solution to the linear Schrodinger equation, 
but in the effective potential V e ff = V — bu 2 . The other way boils down to picking u that ensures that 
(27) leads to a normalizable and positive R 2 . Inserting now R and u into (28) gives us the potential 
that supports this configuration. Nevertheless, the solutions found in these ways may not necessarily 
yield to a simple physical interpretation. For instance, by employing the second method for c = 0, 
a = —2, and u(x) = 1 + x 2 one obtains R 2 = 1/(1 + x 2 ). This configuration is obviously well localized 
and thus normalizable. It exists in the effective potential V e ff = (2x 2 — 1) /2(x 2 + l) 2 — 7 (1 + x 2 ) 2 
which is unbounded from below as 7 is positive. The energy Eqm of this stationary solution is finite 
and equal zero even if the potential V e jj can be negative in the whole domain for a sufficiently large 
7. However, its field-theoretical energy is infinite. The configuration in question does not affect the 
standard Ehrenfest relations. 

Similar solutions appear in the linear equation as well. For example, one can obtain exactly 
the same solution as considered above in the linear case if c 7^ 0. In the modification discussed, 
the solutions in question might, however, play more important and perhaps a novel role due to a 
more pronounced contribution of the phase to the equations of motion. The major problem with 
these solutions in both linear and nonlinear theory is that they seem to be vulnerable to generic 
small perturbations. The perturbed configurations can have infinite energy and therefore cannot be 
considered physical. This implies that the functional domain of stability of the solutions discussed is 
very small. In fact, they may not be stable at all, that is, an infinitely small perturbation can destroy 
them. Such a statement would require a formal proof, but since we cannot provide it here, we offer 
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it as a conjecturej] Consequently, it is safe to say that no other normalizable stationary solutions of 
physical interest beyond those to the linear equation appear in the Staruszkiewicz modification. In 
principle, one can also address the issue under consideration as a spectral problem of equation (9). 
However, in the case of nonlinear equations, each spectral problem requires a separate mathematical 
analysis as, unlike for linear operators on Hilbert spaces, there is not any general theory that could 
tackle such problems. 

One of the most general frameworks in which to seek solutions to the modified Schrodinger equation 
is that of dynamical systems. For simplicity, we will employ it here only for the one-dimensional case, 
in which one has 

, 1 



uj = — 

7 



<9e 1 fn . N 
— + -e (2( V + v) 
at m 



(' = 2m 



'dS r „ , 1 2 
- + V (x,t) + - V 2 



-e 



(29) 



(30) 



2eC, 



V 



v 



v. 



(31) 
(32) 
(33) 

with the variables rj, e, (, u, and to defined as r\ = S' = dS/dx, e = R 2 , ( = R' / R, v = S", and uj = S"'. 
It should be noted that the dynamical system (29-33) "evolves" in x, the actual time variable t being 
a parameter. 

Looking for stationary points of this system, we dismiss e = as unphysical, which leaves us with 
an alternative ( = uj = v = and de/dt = 0, and dS/dt + V(x,t) + rf /2m = 0. As can easily be 
verified, the plane wave described by S = —Et + kx and R = 1, with the normalization condition 
realized by imposing the periodic boundary conditions, is a stationary point of this dynamical system.^ 
In fact, it is probably the only such a point. 

As a way to demonstrate the capacity of this method let us now find out how the other solutions 
to the nonlinear Schrodinger equation discussed in the present paper can be identified within this 
framework. The coherent wave packet with R 2 = e and S given by (20-21) corresponds to 77 = 
— \/2|o;| sin(c<jt — 8)/xq, ( = — (x — Xq\/2 cos(iot — 5)^j /xq, and v = u = 0. One deals here with an 
effectively two-dimensional phase space (for any fixed t) spanned by ( and e as they are the only 
variables allowed to change in "time" x. For comparison, the Gaussian wave packet with R and S 
specified by 



R 



mto 



7T (f + t 2 ) 



1/4 



exp 



mtox 



2 (t 2 + t 2 ) 



and 



S 



mtx 



2 it 2 + tl) 



1 t 
arctan — . 

2 t Q 



(34) 



(35) 



lives in an effectively three-dimensional phase space. Indeed, in this case the motion takes place in vari- 
ables r] = mt\x\/ (t 2 + tl),( = —mt \x\/ (t 2 + t 2 ^, and e, the other variables being v = mtj (t 2 + tl) 

7 Since the solutions in question appear in both linear and nonlinear theory, we will discuss them in greater detail in 
a separate paper hoping to resolve this issue in it [ p^[ . 

8 This observation was first made by A. Z. Gorski and P. O. Mazur in 1984 |17| . 
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and u) — 0. They both are constants for the actual time fixed. With the time allowed to run, the man- 
ifolds on which the evolutions of these packets occur are three- and four-dimensional, correspondingly. 
This comparison of ordinary wave packets with the coherent states is just another way of demonstrat- 
ing that the phase space of the latter is smaller than that of the former. It is due to this very fact 
that the coherent wave packets are also called "the minimal phase space packets". The KdV soliton 
does not represent any particular plane in this phase space. 



5 Conclusions 

The paper presented was aimed at recalling and a further discussion of the nonlinear modification 
of the Schrodinger equation proposed by Staruszkiewicz. Historically motivated by a longitudinal 
modification of the electromagnetic action, here it has been treated as postulated on the grounds 
of physical consistency and certain attractiveness similarly as other notable nonlinear modifications 



of this equation [T0|, [18| . Unlike these, the Staruszkiewicz modification does not introduce any new 
dimensional constants in the system of natural units as long as it is formulated in three dimensions. 
As argued, due to the energy ambiguity typical of nonhomogeneous modifications of the Schrodinger 
equation, depending on which physically relevant properties one decides to maintain, one can come 
up with three different formulations of the Staruszkiewicz modification, although we believe that 
the arguments presented favor the moderate formulation which assumes that the correct definition 
of energy is given by the expectation value of the Hamiltonian and rejects any other definitions as 
spurious. It is within this formulation that the largest number of solutions can be found, including 
the Gaussian wave packet that could not be incorporated in the original version of the modification. 
The formulation in question admits also some other physically interesting solutions, among them the 
coherent state for the potential of harmonic oscillator, shared though by all the formulations, and the 
KdV soliton for another suitably chosen but time-dependent potential. None of these solutions alters 
the standard Ehrenfest relations. 

One can speculate that Staruszkiewicz's modification describes some part of reality that due to 
rather an elusive nature of the phase of wave function has escaped our realization. In fact, until 
relatively recently the phase had been treated as an object of rather secondary importance to the 
understanding of quantum-mechanical systems. However, it turned out that in a number of physical 
situations such a position could not be maintained for it would lead to an incomplete physical de- 
scription. The phase of the wave function encodes relevant information on the quantum evolution as 



manifested in the Aharonov-Bohm effect |T3| and other phenomena of similar nature [IS, p0[ (see also 



pT| for a more complete list of references). Another motivation to study this modification stems from 
the long-standing problem of the collapse of wave function. Since it is rather reasonable to expect 
that the phase plays an important part in this process, one hopes that the modification discussed can 
offer some insight into the physics of this phenomenon. 

In all its formulations, the Staruszkiewicz modification introduces very minimal changes to linear 
quantum theory. Not only are the stationary states and thus also the atomic structure unchanged, but 
also the energy of nonstationary phenomena remains the same in most if not all physically relevant 
cases. Despite these, the modification in question differs from linear quantum mechanics in its classical 
limit. It is, however, not out of the question that to solve the riddle of the measurement problem and 
perhaps to understand other phenomena currently viewed as anomalous, the idea of modifying the 
classical limit of quantum theory must be seriously considered. Being a fundamental theory of nature, 
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it is the quantum theory that should tell us what the world is like in the classical limit. Imposing 
such a limit on the basis of our apparently "classical" collective experience may deprive us of the 
understanding of phenomena that are not so common to this experience and thus, sometimes, not 
even widely regarded as real. As noted by Einstein, it is theory that is supposed to tell us what is 
observable. The Staruszkiewicz modification offers the subtlest conceivable proposal to extend the 
classical limit. Moreover, as partially indicated by this non-standard limit, it is very plausible that this 
modification is not linearizable and therefore it could describe phenomena that cannot be captured 
by the linear theory. 
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